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Abstract 

A rotating bosonic many-body system in a harmonic trap is studied with the 3D-Cranked 
Hartree-Fock-Bogohubov method at zero temperature, which has been apphed to nuclear many- 
body systems at high spin. This method is a variational method extended from the Hartree-Fock 
theory, which can treat the pairing correlations in a self-consistent manner. An advantage of this 
method is that a finite-range interaction between constituent particles can be used in the calcula- 
tion, unlike the original Gross-Pitaevskii approach. To demonstrate the validity of our method, we 
present a calculation for a toy model, that is, a rotating system of ten bosonic particles interacting 
through the repulsive quadrupole-quadrupole interaction in a harmonic trap. It is found that the 
yrast states, the lowest-energy states for the given total angular momentum, does not correspond 
to the Bose-Einstein condensate, except a few special cases. One of such cases is a vortex state, 
which appears when the total angular momentum L is twice the particle number N (i.e., L = 2N). 
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I. INTRODUCTION 



Since the Bose-Einstein condensation (BEC) was realized in trapped dilute atomic gases 
at ultra-low temperature, theoretical studies of the BEC have been rapidly developing. 

In the early stage of study, ultra-cold alkali atoms such as ®^Rb and ^^Na were mainly 
used for a formation of the BEC. Many theoretical analyses were performed with the Gross- 
Pitaevskii (CP) equation, in which the two-body interaction takes a delta-function form. 
In fact, the two-dimensional CP equation is suitable particularly for the study of the BEC 
made of alkali atoms in a cylindrical trap, where the s-wave scattering is dominant. 

Thanks to recent experimental developments, non-alkali atoms and molecules were also 
cooled down to form the BEC. Such an example is seen in the condensate of Cr atoms [1]. The 
inter-atomic potential in the BEC made of Cr atoms cannot be approximated exclusively 
by the delta function because of the strong dipole moment carried by a Cr atom j^, J 
Therefore, the dipole-dipole interaction (including a tensor force) could be also responsible 
for the many-body dynamics, which involves the d-wave scattering in the BEC. Another 
example is the condensate of molecules |4|, where the delta function is not appropriate for 
the intermolecular interaction because of the anisotropic nature of the interaction. 

The BEC can be rotated with the state-of-the-art experimental techniques, such as the 
laser spoon js^ . Owing to these techniques, it was demonstrated that the BEC undergoes the 
quantum phase transition to vortex states. The "cranked" CP equation was applied to the 
analysis of the rotating BEC, and it successfully explained the quantum phase transition, 
including the triangular lattice of vortices [sl, 3]- In the early stage of the study of the 
rotating BEC, the two-dimensional CP equation was mainly used for theoretical analyses 
8j, because the BEC was only rotated about the fixed axis. 

Recently observed phenomena, such as precession 9, 10] and bending [ll] of vortices, 
require the rotational axis to move around in a time-dependent manner, with respect to a 
certain coordinate frame. These phenomena attract much interest in terms of the three- 
dimensional spatial structures of the vortex states. Hence, the three-dimensional CP equa- 



tion is now being applied to various vortex states 
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In addition, a topo- 



logical technique was deve 
magnetic field of the trap 



oped to produce vortices with spin 2 or 4 (h), by reversing the 
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18| . (Below, we take the unit for angular momentum to be 



h = 1.) The formation mechanism of such vortices requires three-dimensional motion of the 
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vortices. In these situations, the total angular momentum vector needs to be treated in a 
three-dimensional manner. 

In this way, everytime new experimental progresses are achieved, new physical situations 
are created, to which the original GP equation cannot be applied in a naive way. The 
cranked Hartree-Fock-Bogoliubov (cranked HFB) method, which has been used to describe 
rotational states of atomic nuclei Q, could be a useful and powerful approach in order 
to deal with these new situations. 

Rotating BEG systems in a trap has been analyzed with the GP equation, or the exact 
diagonalization method using a set of the truncated basis. Mottelson was the first to discuss 
the "yrast" structure [21 1. Here, the yrast states mean the lowest-energy states for given 
angular momentum. He proposed a scenario 22| that the quadrupole excitation is dominant 
when the total angular momentum L is much less than the particle number N {L << N), 
and that all the bosonic particles will occupy the p state when the total angular momentum 
becomes equal to the particle number. This situation can be interpreted as a creation of 
a vortex state. This prediction was numerically verified by himself and his collaborators 
using the two-dimensional GP equation Bertsch and Papenbrock also verified this 



model Hamiltonian 
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26|. 



27|. 



prediction using the diagonalization of the two-dimensiona _ 
Further detailed analysis was performed by Nakajima et. al. 

In this paper, we apply the cranked HFB method to a simple schematic model, where 
bosonic particles interact weakly through the repulsive quadrupole-quadrupole interaction. 
This model is too simple to describe the detailed structure of realistic systems. But, in 
limited situations, the quadrupole-quadrupole interaction becomes a phenomenologically 
valid interaction that can reflect physically essential properties of the realistic interaction. 
For example, our model can describe low-energy rotational excitations of weakly interacting 
alkali atoms, as discussed by Mottelson 22]. Then, using the density matrix of the yrast 
states, as well as its eigenvalues and eigenstates, we compare our results with the other 
methods. Greation of a vortex state is also discussed within the framework of our model. 

In Sections II and III, we present how the cranked HFB theory is extended so as to 
calculate not only fermionic systems but also bosonic ones. Unlike the GP equation, we do 
not assume the inter-atomic potential to be the delta function. Also, we do not suppose an 
a priori existence of the condensate, which is the essential assumption in the GP equation. 
The cranked HFB theory is a constrained mean-field theory, and the value and direction of 
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the total angular momentum vector are controlled in the calculation. With this method, it is 
expected that we can numerically analyze not only structure of dilute many-boson systems 
in a trap, but also superfluidity produced by ultra-cold many-fermion systems in a trap. 
In the present study, we focus on the study of weakly interacting Bose systems, and an 
application of the cranked HFB theory to a simple model is presented in Sections IV and V. 



II. CRANKED HARTREE-FOCK-BOGOLIUBOV THEORY 



We describe the cranked HFB method that can be applied to rotating particles interacting 
two-body interactions. This method was originally proposed for a description of nuclear 
rotation 



19 



20[ |. but we extend the method to deal with not only fermions but bosons. 



Let cJj,,Cq, be the creation and annihilation operators of the single particle state {^{a) = 
■^0(0 5 where ^ represents the real-space coordinates, spin coordinates and nuclear spin of 
particles. The creation and annihilation operators of quasi-particle al,ai is given by the 
Bogoliubov-Valatin transformation 28l. l29l|. 



a 

O'i = ^ ] U^j^Ca + Vai'^a 



(1) 

(2) 



The operators a^, Ca obey the following commutation rule. 



[ai,a|]± = 5ij, [a\,a\]± = [ai,aj\-_ 



(3) 



± 



0, 



where the upper sign (-I-) applies to fermions and the lower (— ) to bosons. To satisfy the 
commutation rule ([3]), we need the following relations. 



U^U ± 1/V = 1, U^U* ± V^V* = 1, 
[/V* ± V^U* = 0, V^U ± U^V = 0. 



(4) 
(5) 



Based on the variational principle, the U and V are determined. The variational ansatz is 
chosen to be 

m=Njexp{J2bapcicl)\0), faP = J2^a^{U*-'),p, (6) 



a/3 



4 



where, Nf is a normalization constant and |0) is the true vacuum. The variational state |<l>) 
corresponds to the vacuum in the quasi-particle basis, that is, = 0. The many-body 

Hamiltonian including a two-body interaction V{C,i,C,2) is generally written as 

where Va/S'yS is given by 



The one-body part H^^ includes the kinetic energy and the spherical confinement potential. 
Using Wick's theorem, we can represent the Hamiltonian by the total energy E and the 
one-body Hamiltonian h, as 



H = E + h + ^^ Vaf3^s : clc^csc^ :, (9) 

afS-fS 

E = + + (CI4)(C5C^) + 2(4c^)(4c5)), (10) 

h = ^H^p- clcp ■ +i ^ V«/3^5((C5C^) : clcl : +{clcl) : c^c^ : +2(cJ,c^) : c^Cs :),(11) 

where : ■ ■ ■ : is the normal order product with respect to a, and an abbreviated expression 
is introduced for an expectation value, (O) = 

The HFB wavefunction |$) is determined through the variational principle with con- 
straints Cn, 

5mH-Y,f^nCnm=0, (12) 

n 

where /i„ is a Lagrange multiplier. In this method, three components of the total angular 
momentum and the particle number are constrained. Further constraints are imposed on 
the following quadrupole operators, Bi = \J^yz, B2 = sj^^^^ = \l^xy, in order to 
fix the intrinsic coordinate axes of the system along the principal axes of the quadrupole 
moments. Therefore, we have seven constraints in our calculations. 

($141$) = J,, ($1 J,|$) = J„ ($1 J.|$) = J., (13) 
($|Bi|$) = 0, ($|52|$) = 0, ($|53|$) = 0, (14) 
($|iV|$) = N, (15) 



In Eq.(IT2l). these constraints are represented with C*„'s as Ci = Jx,C2 = Jy^C^ = Jz,C4^ = 
-Bi,C*5 = -82,(76 = B^jCj = N. In particular, the term —Yln=if^n^n been called the 
"cranking term" in nuclear high-spin physics, because it simulates the effect of the Coriolis 
force in the rotating mean-field system. 



III. METHOD OF STEEPEST DESCENT 



As mentioned earlier, we determine the HFB states, following the variational principle. 
By multiplying a unitary operator to an arbitrary initial HFB state |$), another HFB state 
1$') is obtained. This transformation is considered as a variational procedure with respect to 
the matrices U and V of the Bogoliubov-Valatin transformation, Eq.(l2]). The transformation 
is iterated until a local minimum is found to satisfy Eq. ffT^ . 

Now, let us explain how the unitar y op erator is given within our framework. First of all. 
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3l|, the unitary transformation of the HFB state 



from the extended Thouless theorem 
is expressed as 

|$')=exp(f/)|<|.), (16) 
where d is an anti-Hermitian operator d = —d\ which is generally expressed as, 

A quasi-particle basis is then transformed in the following way. 

exp(±P)^ I , (18) 



a'} \ / e-^a]e^ ^ ' 




where 

D = I " " I , (19) 

yd* J 

and d is a matrix representation of d. We choose the anti-Hermite operator d as, 

d = [rir + s]^ (20) 
where we define [O]'^ = ^[Af,0] and the quasi-particle number operator is given as A/" = 



6 



J2i operator s and the single-particle Routhian f are respectively defined as 

S = Y.^nC^^ (21) 

n 

f = /i-^/i„C'„. (22) 

n 

Then, these parameters rj, 5n, and /i„ are determined through a minimization of 
under the constraints (^ICjl^) = q. The parameters 6n and /i„ are evaluated by expanding 
($'|C*j|<I>') up to the first order in 5„ and fin- That is, 

h = $^i.,t(c.-(<f|Q|$)), (23) 

i 

where L^i = {^\[Ci, [C*fc]"]|$). Whereas, the parameter i] is determined from a minimization 
condition for E = {^'\H\^'), through expanding E up to the second order in rj. As a 



consequence, we have 



($i[[^,[f]'^],[f]'^]i<i>) 



To check the convergence for the self-consistency in the calculation, it is convenient to 



define the norm of d as \d\ = y iMijl^- our calculations, a criterion for the convergence 
is given whether \d\ is less or greater than e = 1.0 x 10~^. When \d\ < e, we judge that 
the convergence is numerically achieved. Otherwise, the iteration for the self-consistency 
continue until \d\ meets the above condition. 

IV. A SCHEMATIC MODEL 

To examine the convergence procedure of our method (3D-cranked HFB) for the bosonic 
case, let us consider a toy model, which is similar to a realistic system of dilute ultra-cold 
Bose gases confined by an isotropic harmonic oscillator potential V{r) = ^Muj'^r'^, where 
M represents the atomic mass. We choose the single-particle states to be the harmonic 
oscillator states, {r\a) = (r|cj^|0) = ('")«'°^«QmQ (^0); that is, a product of the Laguerre 
polynomial and the spherical harmonics. Let us denote this basis as a = {ua, la, ^a)- 

Mottelson discussed that the quadrupole correlation is important in the low angular 



momentum region 



22| . In accordance with his proposition, the following Hamiltonian is 



considered in our calculations. 

1 ^ 

i^modci = Ho + -kY, i-rQ-f.Qf. + gP^P, (26) 

Ho = J] (2^- + ^° + 2)^'^^-^«' (27) 

a 

Q, = J2Wr'ClP\P)cic^, (28) 

P = A/2/a + IjlgmJa - ma\00)CaCa, (29) 

a 

where Ca is the annihilation operator corresponding to the state a = {na,la, —i^a), and 
Ci'^\Q) is related to the spherical harmonics through Ci''\fl) = sj^^^^kni^)- This Hamil- 
tonian should be considered as a simple model for weakly interacting dilute atomic gases. 
The parameter k represents strength of the quadrupole-quadrupole interaction, and k is 
positive in this work, to treat the repulsive two-body interaction. The last term of Eq. (!26ll 



is called the pairing interaction. [24] The parameter g represents the strength of the pairing 
interaction. In the present calculation, we set g to be zero, for the sake of simplicity. The 
one-body Hamiltonian (fTTl) in this model is given as, 

^model = ^0 + ^ $^(-)^(Q-,.)Qm + 9{P^)P- (30) 

The oscillator energy for the isotropic harmonic oscillator states, hw, is set to 1.9 (meV). 
The quadrupole-quadrupole interaction k/A^ is repulsive and its strength is set to 0.1, 0.5 
and 1.0 (meV//im^), where A is the mass number of an atom. The single-particle model 
space used in this calculation are the harmonic oscillator states of the Os, Is, Od, 2s, Id, Og, 
Op, Ip, Of, 2p, If and Oh states. 

To prepare the initial state, we use the deformed quadrupole mean field. 



Wm = i^o - -M^V^ (^/?cos7C'(^) + ^/^sin7(Cr (r]) + Ci^2H^))J • (31) 

The quadrupole parameters (/3, 7) are useful measures to think about the shape of the 
many-body system. /3 is a measure for elongation or stretching, while 7 for triaxiality or 
deviation from axial symmetry. For example, a spherical shape has /? = and nuclear 
superdeformation has typically /3 ~ 0.6. The triaxial parameter 7 gives axial shapes when 
7 = 0° and 7 = 60°. The former shape corresponds to the so-called "prolate" shape, which 
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is similar to a kiwifruit, while the latter to "oblate" shape, similar to a mandarin orange. 
By definition, a triaxial shape is invariant with respect to an operation: 7 — > 7 + 120°. 

We diagonalize the Hamiltonian (I3T1) for (/5,7) = (0.01,0) and obtain the single-particle 
states. The initial state is created as a Hartree-Fock state, that is, all the single-particle 
levels are occupied below the Fermi level. 

Since the initial state is symmetric under rotation about the z-axis, collective rotation 
about the z-axis is suppressed. In other words, we cannot crank the state around the symme- 
try axis. To induce angular momentum, we initially set the constraints of the total angular 
momentum to {Jx,Jy,Jz) = (0.05,0,0). As a next step, we tilt the angular momentum 
vector to {Jx, Jy, Jz) = (0, 0, 0.1). With this procedure, we can increase Jz up to 20, with a 
step, AJz = 0.05. 



V. RESULTS AND DISCUSSIONS 

Figure [H represents the total energy E as a. function of the total angular momentum L. 
Four lines are plotted in this figure. One dotted line represents huj^L + ^N) and the other 
three lines are calculated with the different quadrupole-quadrupole interaction strengths, 
which are k/A^ = 0.1, 0.5 and 1.0 {meV / fim'^) . Despite the different values for n, These 
lines are nearly identical. In other words, E is almost independent of k. 

We also find that E increases almost in proportion to L, that is, E oc L. This result can be 
explained from a microscopic point of view. As L is increased, single-particle excitations are 
induced one by one through the two-body interaction, so as to satisfy the angular momentum 
constraints. In other words, AE = huAL, where hw is the single-particle energy spacing of 
the isotropic harmonic oscillator. Due to the quantum statistics for bosons, this excitation 
mode can continue until the number of particles occupying the ground state becomes zero. 



This linear behavior is already noticed by other authors 25|. However, with a careful look 



at our numerical result, there is a slight deviation from the linearity in the the total energy 
(hhj{L + |A^)). This small deviation is caused by deformation of the mean field, which has 
a role to mix the single-particle orbits. The evolution of the quadrupole deformation in 
response to rotation is discussed below. 

Figure [2] shows the deformation parameters j3 and 7. The strengths of the interaction are 
set in the same manner as in Figure 1. These deformation parameters are self-consistently 
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calculated as 



7 = tan 



1 



(Q22) 



(32) 



(33) 



V2{Q2o)' 



The unit for the quadrupole moment {Muj^/k) is derived from the consistency at L = 
between the one-body Hamiltonian (l30il and the deformed mean- field Hamiltonian (13T1) . The 
figure shows that the deformation parameters do not depend on the interaction strength very 
much, although 7 shows minor differences at low spin. When the total angular momentum 
is small {L < 1), the mean field has a almost spherical shape. This is because the trapping 
potential is spherical and the present two-body interaction is repulsive. As L is increased, f3 
increases gradually. In the small angular momentum region (L < 5 ), 7 is not 180°(= 60°). 
(See the right panel of Figure O) This result means that the shape of the mean field is 
not axial-symmetric, but triaxial. When the quadrupole-quadrupole interaction becomes 
stronger, the deformation tends to prefer a more triaxial-deformed shape. However, any 
of the three cases ends up with the oblate shape (7 = 180°) at high angular momentum 
(L > 5). A reason for this tendency can be explained as the following: When the quadrupole- 
quadrupole interaction is strong, the harmonic oscillator states having the different magnetic 
quantum numbers become more mixed through the interaction. As a result, the magnetic 
quantum number is no longer a good quantum number. This is nothing but axial symmetry 
breaking, or an emergence of triaxiality. It should be noted, however, that 7 is substantial 
only at low angular momentum, where (3 is very small. In other words, when the elongation 
is small (/5 ~ 0), the triaxial degree of freedom is irrelevant in terms of a deviation from a 
spherical shape. That is, in our calculation, the shape of the mean field can be regarded 
to be almost spherical in the small L region. On the other hand, for the higher angular 
momentum region (L > 5), 7 is almost constant to be 180°, meaning that the mean field 
becomes an oblate shape. 

The condition for the BEG of weakly interacting bosonic atoms in a trap is given by 
Nv/tijjj « 1 [22], where v is an expectation value of the two-body interaction, while 
huj represents the single-particle level spacing. In our calculation, Nv corresponds to an 
expectation value of the quadrupole-quadrupole force, that is, /t(((52o)^ + 2(Q22)^)- The ratio 
Nv/huj is then estimated to be 4.4 x 10~^(3A'^ / k. According to our calculation, deformation 
is up to /5 < 1, so that the ratio is of order of 10~^ to 10^^ for our three choices of the 
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interaction strength {n/A^ = 0.1, 0.5, 1 meV/fim^). This result means that our calculations 
can be regarded as a weakly interacting many-boson system. 

Figure [3] shows the occupation probability, paa for k/A^ = 0.1 meV/fim^, where pa/s is 
the density matrix defined as 



{V is a matrix appearing in the Bogoliubov-Valatin transformation, Eq.([2]).) Although 
the occupation probabilities for 1 meV/fim^ are not plotted in Figure El we have 

calculated these occupation probabilities and found that they are almost same as that of 
k/A^ = 0.1 meV / fiTo^ . This result indicates that the wave-function does not strongly depend 
on the strength of the quadrupole-quadrupole interaction. 

At low L, the (OsO) states are the major component in the HFB state. The higher the 
total angular momentum, the more the {0d2) state admixes with the (OsO) state. The {0g2) 
component is also mixed at L ~ 2N = 20, while the (OsO) component vanishes. This result 
suggests that the yrast state changes its structure gradually. 

Figure m shows the eigenvalue Ua of the density matrix paf^ , as a function of L. Only the 
plot for k/A'^ = 0.1 (meV//im^) is displayed because the occupation is almost independent 
of k/A"^. The largest eigenvalue in the figure is equal to the total particle number N = 10, 
and this situation happens only at L = and 2N{= 20) In these cases, all the particles 
occupy only one single-particle state, which can be regarded as the condensate state. On 
the other hand, between L = and 2A^ (0 < L < 2A^ = 20), the particles are shared by the 
two states and ip^, where ua + i^b = N = 10. This result indicates that most of the yrast 
states are non-condensates, but a mixture of two single-particle components, t/'^ and ip^. 

In Figure [5l the eigenstates and ip^ are decomposed into the single-particle basis, 
and their components are displayed in terms of probability (f^)^, where {ip"") = Xlo^al'^)- 
First, as shown in the right panel of the figure, the state at L = is found to have a 
condensate structure into the (OsO) state, which is consistent with the occupation number 
calculation shown in Figure [3l As L is increased, the (OdO) component starts to mix with 
the s component although the contribution from the d state is minor. This mixture is 
consistent with the growth of the quadrupole deformation, as shown in the left panel of 
Figure [21 When L < 5, axial symmetry is broken around the cranking shown 
in Figure [21 In this situation, the whole system rotates in a collective manner, which is 




(34) 
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consistent with Mottelson's model claiming that the yrast structure is dominated by the 
collective quadrupole excitation at low L. Such a collective rotation was actually observed 
experimentally at lower rotational frequency before vortices are formed However, the 
linear dependence of on L at low L seen in Figure [1] implies that the collective mode is 
not the major mode in our model, but that the single-particle excitations are. Next, in the 
right panel of Figure 0, the state iIj^ is decomposed into the d and g states. This is because 
the quadrupole-quadrupole interaction does not mix the states having different parity, so 
that the particles in the s state can not be excited to the p state, but the d or g states. The 
major component is the (Oc/2) state when angular momentum is small {L < 10), whereas the 
{0g2) state starts to mix in the high angular momentum region (L > 10), due to the onset 
of deformation in the mean field. Although the (0(74) state is included already at L ~ 0, this 
component can be considered as a minor component because the occupation ub is nearly 
in this region. 

According to the previous discussion, at L = 2A^, the many-body state goes into the 
condensate in the state ip^, which is a mixture of the {0d2) and {0g2) states. Both of these 
{0d2) and {0g2) states have angular dependence sin^ 9, so that the density along the 2;-axis 
vanishes for these states. Therefore, the condensate at L = 2A^ is considered as a quantized 
vortex state carrying 2 (ti). 

The HFB solutions obtained for the intermediate L values (0 < L < 20) are quite different 
from the solutions obtained with the GP equation. This is because the GP equation assumes 
an a priori existence of the condensate for the whole range of L, which can be expressed 
as a Hartree state, (aQ)^|0). Although our HFB ansatz includes this condensate state as 
a special case, the mathematical form for the HFB state is generally more complicated to 
allow a linear combination form of multiple Hartree states, in accordance with Eq. iQ. If 
the condensate is realized at any L, all the particles should occupy the one single-particle 
state that is generally expressed by a liner combination of the basis states, such as the s, d 
and g states. 

As shown in Figure HJ the yrast structure changes smoothly from the ■i/'^-dominant states 
to the '?/'^-dominant state, as L is increased from to 2A^. Considering that the latter state 
at L = 2A^ corresponds to a vortex state, a formation of the vortex starts as a shallow dent 
in the center at small L. As the amplitude of ip^ becomes larger for increasing L, the depth 
of the dent becomes deeper. Finally at L = 2A^, a complete vortex is formed, where density 
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becomes zero. This mechanism of the vortex formation is very different from the one derived 
from the GP equation. In the calculations using the GP equation Q], a vortex enters from 
the "outside" of the system due to a continuity of the many-body wave function. This sort 
of process needs a odd-number multipolarities, such as a dipole (A = 1) and an octupole 
(A = 3) correlations, which are missing from our present model. 

These higher multipole correlations may play an important role also in a formation of 
vortex lattices, which violates rotational symmetry of the system. In the present framework, 
vortices appear always in the center of the system due to the axial symmetry possessed by 
the system. To allow multiple vortices to appear away from the center, we need, at least, a 
mixing between these states of (OsO), {0d2) and {0g2) to break the symmetry. This situation 
is realized only when {Q2±2), which is contained in the single-particle Hamiltonian, Eq.f l30l) . 
is non-zero. In other words, the deformation parameter 7 should not be equal to 0° or 180° 
to allow triaxial deformation bringing anisotropy to the system. However, in the present 
calculations, 7 takes the value of 180° in a wide range of angular momentum (5 < L < 20). 
Consequently, the vortex lattice is not produced in our model. 

VI. SUMMARY 

Extending the 3D-cranked Hartree-Fock-Bogoliubov method, we have performed the nu- 
merical calculation for a rotating many-boson system interacting through a weak and re- 
pulsive interaction, trapped inside an isotropic harmonic oscillator potential. Unlike the 
Gross-Pitaevskii equation, our calculation does not assume the existence of the condensate 
a priori, but general many-body states in the framework of the HFB method. 

We applied the method to a simple model where the two-body interaction is chosen to be 
a separable type called the quadrupole-quadrupole interaction. Parity is conserved for this 
interaction, so that only Al = 2 excitations is allowed through the quadrupole-quadrupole 
interaction. 

First of all, at L = 0, our calculation shows that the HFB state is interpreted as a 
condensate into (OsO). 

As increasing the total angular momentum L from to 20, the particles in the yrast state 
transfer from (p^ to (p^. which have the magnetic quantum numbers = and = 2, 
respectively. At low L, triaxial deformation is formed to allow collective rotation around the 
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cranking axis. However, the linearity E cc L implies that the major excitation mode is still 
single-particle excitations. At higher L, the system becomes oblate, that is, axial symmetric 
around the rotating axis. Angular momentum is thus produced by migrating from s state 
to d state. These single-particle excitations are induced from the s to the d state through 
the quadrupole-quadrupole interaction, in agreement with the prediction by Mottelson. 

Finally, at L = 2N, the HFB state becomes another condensate, in which all the particles 
occupy the single state expressed by a linear combination of the (0d2) and {0g2) states. This 
state can be interpreted as a vortex state having angular momentum 2{h). 

In this way, the yrast structure changes gradually and smoothly in our framework. This 
result can be accounted by the following two effects: One is our choice of the HFB ansatz 
which allows a linear combination form of multiple Hartree states. The other is the finite- 
number effect of the total particle number {N — 10). These effects surely needs further 
investigations in the future studies. 

We are currently proceeding to extend our programming code to deal with more realistic 
inter-atomic potentials, and plan to examine the effect of the pairing interaction in the HFB 
framework. 
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FIG. 1: The total energy (meV) as a function of the total angular momentum for k/A"^ = 0.1,0.5 
and 1.0 (meV//um^). 
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FIG. 2: The left and right panels represent the deformation parameters (3 (unit-less) and 7 (degree) 
as a function of the total angular momentum, respectively. 
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FIG. 3: (Color online) The occupation number p^a of the wave-function for n/A^ = 0.1 (meV/ /xm^). 
The numbers in the figure {nalania) stand for the principal, orbital and magnetic quantum numbers 
of the harmonic oscillator state. 
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FIG. 4: (Color online) The eigenvalue Va of the density matrix p^p for k/A^ = 0.1 (meV/;um^) as a 
function of the total angular momentum. Only the three eigenvalues uq are ploted because 

the other eigenvalues are almost zero. 
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FIG. 5: (Color online) Left (right) panel represents the components for eigenstate ■i/'^ (V'^) of the 
density matrix as a function of the total angular momentum, respectively. The eigenstate V'^ (V*^) 
corresponds to the eigenvalue i^a (z^b) in Figure HI 
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